Symmetry-induced gap opening in graphene superlattices 
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We study nxn honeycomb superlattices of defects in graphene. The considered defects are missing p z orbitals 
and can be realized by either introducing C atom vacancies or chemically binding simple atomic species at the 
given sites. Using symmetry arguments we show how it is possible to open a gap when n = 3m + 1 , 3m + 2 (m 
integer), and estimate its value to have an approximate square-root dependence on the defect concentration x = 
l/n 2 . Tight-binding calculations confirm these findings and show that the induced-gaps can be quite large, e.g. 
~ 100 meV for x ~ 10~ 3 . Gradient-corrected density functional theory calculations on a number of superlattices 
made by H atoms adsorbed on graphene are in good agreement with tight-binding results, thereby suggesting 
that the proposed structures may be used in practice to open a gap in graphene. 



The recent fabrication of graphene 01, a one atom-thick 
layer of carbon atoms arranged in a honeycomb lattice, 
has triggered a wealth of studies in both fundamental and 
applied science. Graphene is a a zero-gap semiconduc- 
tor with a linear dispersion at the Fermi level in which 
low-energy excitations mimic the behaviour of relativis- 
tic massless fermionsJ2, 0] (see J3, Etl f° r reviews). This 
gives rise to a number of interesting phenomena such as an 
anomalous quantum Hall effect^ |3|] and quasirelativistic 
Klein tunneling [6J. Its unconventional transport properties 
with ballistic transport on submicrometer scaley7|] and with 
carrier mobilities up to 2 10 5 cm 2 V~ 1 s~ l (8|] offer the pos- 
sibility of high-performance interconnects in an hypotheti- 
cal carbon-based nanoelectronics. However, since conduc- 
tivity cannot be turned completely off, pristine graphene 
cannot be used as a transistor in logic applications, where 
high on/off ratios are required|0[T3]. Field-switching ca- 
pabilities depend on the presence (and size) of a gap in 
the electronic structure. Electron confinement can be used 
to open a gap inversely proportional to the confinement 
length when rolling up graphene into single-walled nan- 
otubes or cutting its edges to form nanoribbons. Both these 
possibilities have been exploited, and promising carbon 
nanotube/nanoribbons field-effect transistors realized[9]. 

In this Letter we explore a different possibility made 
available by recent progresses in patterning graphene with 
lithographic techniques! Ill[l2ll . Specifically, we study nxn 
honeycomb superlattices of defects on graphene (defined 
to be periodic structures of defects arranged to form a hon- 
eycomb lattice commensurate with the substrate) and use 
symmetry arguments to design semiconducting structures. 
We show that a gap can be opened by preserving graphene 
symmetry. We estimate that an inverse proportionality to 
the (super)lattice constant n approximately holds for the 
gap size in these structures, and use tight-binding andfirst- 
principles calculations to validate these predictions. The 
computed gaps are sizable and compare favourably with 
those found in nanoribbons at the same length scale [13]. 
In the following we first present our symmetry arguments 
and derive general rules to open a gap in the considered 



superlattices. We then estimate its value and present the 
results of numerical calculations. 

Graphene's unconventional electronic properties are 
strictly related to its Dg/, point symmetry. The A;— group at 
the K (K 1 ) high-symmetry points (D^) allows for doubly 
degenerate irreducible representations, and Bloch func- 
tions built with p z orbitals of the A and B sublattices span 
one of its two-dimensional irreducible representation (ir- 
rep), namely E" . This is enough for the % — %* degeneracy 
at the K(K') point and for the unusual linear dispersion 
relation, irrespective of the level of approximation! 14]. 
These properties are captured by the simple tight-binding 
(TB) model Hamiltonian 

H = -t £ aJbj + H.c.-t' £ ajaj-t' £ b]bj 

<i,j> <ij> «i'J> 

(1) 

where a](b]) is the creation operator for an electron on 
site i of the A (B) sublattice, the first sums run over nearest 
neighboring sites in the honeycomb lattice and the second 
sums run over sites which are nearest neighbors in the tri- 
angular sublattices. The hopping t has been estimated to 
be ^2.7 eV whereas ?' -C t has different values depend- 
ing on the parametrization. When t' = the Hamiltonian 
describes a bipartite system, and we assume that this ap- 
proximately holds for graphene. The consequences of re- 
laxing this approximation will be addressed numerically at 
the end of this Letter. 

Bipartitism has a large impact on the electronic struc- 
ture via the induced electron-hole symmetry. For instance, 
it has long been known that in bipartite systems, at half 
filling, sublattice imbalances due to vacancies strongly af- 
fect the energy spectrum at the Fermi level through the in- 
troduction of midgap states lfl5ll . In graphene, such states 
have a semidelocalized nature with a 1/r dependence on 
the distance from the defect lflrjll . By introducing an equal 
number of defects on each sublattice one restores balance, 
eliminates midgap states and a gap possibly opens. In gen- 
eral, however, there is no guarantee that the gap opens at 
K and does not close somewhere else in the Brillouin zone 
(BZ). Therefore, we focus here on nxn honeycomb pat- 
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terns of defects only, in such a way to constrain (by sym- 
metry) the changes in the band structure and possibly re- 
duce accidental degeneracies. In these structures the high- 
symmetry points where degeneracy is expected are the F 
and K points, that is where the fc-groups (Dg/, and D^i,, re- 
spectively) allow for doubly degenerate irreps. In the fol- 
lowing we use E (A) for a generic two- (one-) dimensional 
irreducible representation, and denote as K„ the K point 
of the nxn superlattice BZ (SBZ). For a strictly bipartite 
system at half-filling, degeneracy at the Fermi level occurs 
when the number of E irreps is odd; when this number 
is even degeneracy, if occurs, can be considered as acci- 
dental. Therefore, introducing symmetric defects in such a 
way to have an even number of E irreps both at the F and 
at the K point a gap generally opens. 

To show that this is indeed possible, we consider a 
generic nxn supercell and count the number of A and E 
irreps generated by the carbon atoms in cell (the so-called 
atomic representations). To this end, it is sufficient to con- 
sider that half of the cell which has D^f, symmetry with re- 
spect to its center / (see Fig JTJ, the remaining half behaving 
similarly. These two half-cells (a and /3 in the following) 
play the role of A and B type of sites in the honeycomb 
superlattice. In each of them, the set of C atoms may be 
grouped in classes of equilateral triangles A, a (A,- a), plus 
a possible atom at I a (Ip) as it happens when n = 3m + 1 
and n = 3m + 2 (m integer). Each triangle spans an A + E 
irrep of the group (centered at 7) which behave as s and 
(p x ,Py) orbitals centered at /; the atom at /, when present, 
spans of course an A irrep. Then, by considering Bloch 
functions built with these s- and (px,Py)~ like orbitals it is 
possible to count the number of A and E irreps for each 
case. At F the Bloch functions built with s-like orbitals 
centered on I a and Ip span two A representations, whereas 
p x ,p y -\\ke functions span two E irreps; at K„ the first gen- 
erate an E irrep whereas the latter span 2A+ E. These 
are also the irreps generated by the p z orbitals of the C 
atoms as long as we discriminate between A and E type 
only. The overall result is given in Table U where the sym- 
bol «3 — O3, 13,23 identifies the three (congruence) classes 
modulo 3, i.e. the sequences n — 3m,3m+ 1 and 3m + 2, 
respectively. 

It follows from TableUthat with the full atomic set (i.e. 
considering pure graphene) degeneracy occurs at the K n 
points when either n £ 1 3 or n £ 23. This is consistent with 
the folding K{K') -> K n (K' n ) and K{K') -> K' n (K n ), respec- 
tively. In the case n £ O3 both K and K' folds to F and there- 
fore a 4-fold degeneracy occurs; this can be considered 
accidental in this context as it cannot be predicted by the 
number of E irreps only. More interestingly, two impor- 
tant results concerning the introduction of p z — vacancies 
are easily proved. 

I. By removing a A{ a , Ajp pair only is not generally 
possible to open a gap. Here, 2(A + E) irreps are re- 
moved both at r and at K n , and no modification occurs 
on the parity of the E sets. Exceptions to this rule are, 
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Table I: Number of irreducible one- (A) and two- (E) dimensional 
representations (per cell) generated by the full atomic basis on a 
nxn honeycomb superlattice at the T and K points of the corre- 
sponding BZ. 



e • • © 




Figure 1: The simplest symmetrically defective superlattice for 
n = 4. Also indicated the two halves of the cell (0!,/3), their cen- 
ters I a ,Ip (red balls) and a triangle Aj a . A defect is a missing p z 
orbital at Ia,Ip- 



of course, those cases where degeneracy is accidental (as 
pure graphene in the O3 case which does show a gap after 
removal of one such pairs, see below). 

II. When n £ I3 or n £ 2.3 removal of the atoms at I a 
and Ip does open a gap. In these cases, the atomic basis 
spans 2A at F and E at K n , thereby turning the number of 
E irreps to be even at both special points. Also in this case, 
exceptions of residual accidental degeneracy are possible. 

The second result provides a very simple way for open- 
ing a gap in two thirds of the cases, that is by introducing 
a p z vacancy at I a and Ip . In practice, this can be realized 
by either removing substrate atoms or using them to cova- 
lently bind simple adsorbates such as H atoms lfl7ll . In the 
latter case, indeed, the C atoms involved in the chemisorp- 
tion process turn their hybridization to sp 3 and effectively 
get out of the n — n* band system. In this Letter we focus 
on the simplest defective superlattices, i.e. those with va- 
cancies at/ a and/^ only, and call them (n,0)-honeycombs. 
In general, the pair of integers (n,p) (p — 0, ..int(n 2 /3)) 
can be used to identify a nxn honeycomb superlattice with 
2p equilateral triangles symmetrically removed from the 
unit supercell, in addition to the atoms at their centers I a 
and Ip when n £ I3IJ23 ■ Fig[T] shows an example, the 
(4,0)-honeycomb. 

It is possible to estimate the size of the induced gap at 
the K n point of these (n,0)-honeycombs. To this end, we 
perform a lattice renormalization[ 18| by making use of the 



3 



K 




Figure 2: Left panel: the renormalized triangular lattice corre- 
sponding to the (4,0)-honeycomb (grey balls). Red and green 
balls mark the position of the A and B defects. Right panel: Bril- 
louin zones for grahene (blue borders) and for the structure shown 
on the left panel (red). Also indicated the four lowest-energy k 
points (K and k, i = 1 — 3) for the calculation in the text (black 
circles). 

bipartite nature of the Hamiltonian H = Hab + Hba- This 
simplifies the problem by halving the state space of inter- 
est. Indeed, since H only allows transitions from the A 
to the B subspaces (Hba) and vice versa (Hba) it is suf- 
ficient to consider the problem in the A space only with 
Hamiltonian Jl9ll Haa = HabHba- For any non-zero eigen- 
value £,■ and eigenvector \YAi) °f this Hamiltonian there 
exist two solutions of the original problem with eigen- 
values £; = ±V5 an d eigenvectors \YA,i) ^ IVb,<)> where 
I Y/fl,;) is denned to be \y/ B .i) = e7 l,2 H BA \ 1//^,); if g< = 0, 
is already a H eigenvector. The converse is also 
true, namely from any eigenvector the two projec- 
tions |VA,i) an d \ Vbj) onto the A and B supspaces satisfy 
Hba \YA,i) = £i\VB,i) and Haa \Ya,i) = ef\VA,i); that is, in 
studying Haa one onrv misses possible zero eigenstates in 
the B subspace[20]. In graphene the renormalized Hamil- 
tonian Haa describes a triangular lattice with on-site en- 
ergy t 2 Z (where t is the hopping term of eq[TJ and Z = 3 
is the coordination number of A atoms in the original hon- 
eycomb lattice) and hopping t 2 between neighbors in the 
triangular lattice. Defects are of two kinds: while A vacan- 
cies translate simply into A vacancies in the renormalized 
lattice, B vacancies modify both the coordination number 
and the hopping between A sites. The renormalized (n, 0)- 
honeycomb for n = 4 is shown in the left panel of Fig|2] 
whereas the right panel of the same figure displays its SBZ 
along with the graphene BZ. The state space at K„ is given 
by its n 2 replicas within BZ and comprises K or K' depend- 
ing on whether n G I3 or n € 23. Quasi-degenerate per- 
turbation theory is necessary to estimate the ground-state 
energy £o(K n ) (and hence the gap E gap = 2^/eq), but it be- 
comes intractable at large n (i.e. at very small defect con- 
centration x = 1 /n 2 ) because a huge number of K„ replicas 
gets close to the K and K' points. Therefore, we consider 
x sufficiently small that the defects are isolated from each 
other, but large enough that a few state calculation is re- 
liable. The smallest set of K„ replicas contains K and the 
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Figure 3: Energy gaps (in units of /) for the (rc,0)-honeycombs 
as functions of x = l/n 2 . Bold black line is the gap at K n as esti- 
mated using Eq. [3] Filled symbols are for results of tight-binding 
calculations with t' = and bold blue lines are best-fit curves (see 
text): squares for n 6 I3, circles for h 6 23; red symbols for the 
gap at K n , blue for the full gap. Open symbols are full gap results 
for t' = 0. 1 1. The inset shows the low-energy band structures of 
the (13,0) (solid lines) and (14,0) (dashed lines) honeycombs as 
obtained from tight-binding calculations with t' = 0. 

k, vectors (2 = 1 — 3) shown in Fig|2] (right panel) for the 
case n—A, and corresponds to the set of Bloch functions 
|i^ ) = \ \if K ) and \ Yt) = |yki) built with p z orbitals at A 
sites; the case n £ 2$ is analogous except that K is replaced 
with K' . I y/o) spans the A' 2 ' irrep of the A:-group at K n (D^), 
whereas {ly 7 ,')} span A'2+E". Thus it is possible to set 
up a two dimensional problem in the A" subspace. The 
corresponding Hamiltonian matrix can be obtained from 
the Haa matrix elements between (graphene) Bloch states, 
Eq|2] 

(wl^AA|%>-f 2 4,k'(3+ J F(k)) 

-xr^k'+g {e ie5A (3 +F(k') +F(k)) (2) 

+^/(k')7(k)} 

Here F(k) = lf =1 e"' k5; , /(k) = lf =1 e"' k5 ' ( where 
8'j and 5, are the vectors joining AA and AB nearest- 
neighbors, respectively ), g is a reciprocal superlattice vec- 
tor and 8a, Sb are the position vectors of the defects in the 
unit supercell. In deriving Eq. |2]periodicity of the super- 
lattice has been used and the defects have been considered 
as isolated (n > 2). With the help of Eq. |2]and of the sym- 
metry properties of F(k) and /(k), the Hamiltonian matrix 
in the above A" space reads as 

A 3 * - xV ^ Fx 1 (3) 

-xV3F x 3+F x - 9x(2 + F x ) 

where F x = F(ki) 3 + 3(27z;/3) 2 x - V3(2^/3) V/ 2 . 

The lowest eigenvalue £ (~ 3xf 2 (0.561 - 0.96 ly/x) ) al- 
lows us to estimate the energy gap E gap at the K n point. 
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Figure 4: DFT results for the energy gap in (n,0)-honeycombs 
made by adsorbed H atoms as functions of n (filled symbols). For 
comparison, also reported are the tight-binding results of Fig. [3] 
(open symbols) with t = 2.7 eV and t' = 0. Squares for n £ I3, 
circles for n e 23. Stars are energy gaps in a number of super- 
lattices with an asymmetric AB pair of defects per cell (ortho- 
and para- dimers not included). The inset shows TB results (nor- 
malized to the number of pairs of defects) for a number of (n, p) 
honeycombs: triangles for p = 1, diamonds for p = 2. (Lines are 
guide to eyes). 



This is plotted in Figj3]as a function of x along with the 
results of tight-binding calculations. The importance of 
including a larger number of K n replicas close to (K, K') is 
evident at small x where the energy gap decreases slightly 
faster than y/x; the best-fit to the numerical results gives an 
exponent ~ 0.66. Tight-binding calculations reveal that the 
minimum gaps occurs either at M or Y (see inset of Fig|3]l, 
depending on x and on the sequence considered, but be- 
have similarly to the gaps at K„. Differences between the 
two sequences appear at large x and reflects the different 
shape of the low-energy bands (not shown). A best-fit of 
£gap = ax a ( 1 + bxP ) to the numerical results, also reported 
in FigEl gives {a,b,a,p) = (3.34,-4.99,0.65,1.14) for 
n G 1 3 and (3.37, +2.9,0.66,0.88) for n S 2 3 . We have 
also investigated the effect of breaking the electron-hole 
symmetry by performing tight-binding calculations with 
t' =^ 0. Introduction of the next-to-nearest neighbors inter- 
action only affects the results at small x, where the valence 
and conduction bands start to overlap at some point x c 
because of the asymmetry introduced in the energy spec- 
trum. The position of this critical value x c shifts to larger 
x when increasing t 1 but remains small for realistic values 
(for t ' = 0. 1 1 used in Fig[3K < 10~ 3 ). Notice that the ef- 
fect of next-to-nearest neighbor interactions is different for 
the two sequences considered. 

To further investigate this point and, more importantly, 
to address the role of electron correlation we performed 
gradient-corrected density-functional-theory (DFT) calcu- 
lations on a number of (n,0) -honeycombs of adsorbed H 
atoms. The details of the calculations are analogous to 
those reported in a previous work[17], except for the fact 



that a finer mesh of k points has been used to correctly 
compute the energy gaps. The results are reported in Fig0] 
up to n = 14 (x ~ 0.005); for larger values of n DFT cal- 
culations become prohibitive. As can be seen from Fig. 
|4]DFT results show a reduced gap size with respect to TB 
ones in all the cases considered but the effect is much more 
pronounced in the 1 3 sequence than in the 23 one. In partic- 
ular, DFT results for n = 5, 8, 1 1 , 14 closely parallel the TB 
ones (f has been set to its accepted value, t — 2.7 eV), de- 
spite the fact that the first refer to a realistic situation where 
defects are H atoms while in TB calculations defects are 
modeled by simple p z vacancies. Discrepancies in the I3 
sequence at small n needs further investigation, though it is 
in line with the introduction of the next-to-nearest neigh- 
bor hopping (Fig. [3J. Similar behaviour was found for the 
gap in armchair graphene nanoribbons where it was ex- 
plained by the modified nearest-neighbor hoppings for the 
sites close to the defect ! 1 311 . In any case, both the size of 
the gap and its dependence on n are promising for future 
applications. 

Finally, to underline the role played by symmetry in 
designing the defective structures, TB results for differ- 
ent superlattices are also reported in Fig. [4] Notice that 
the (n,0) honeycombs considered in this Letter show the 
largest gaps with the minimum number of defects per su- 
percell (see the inset of FigHl where the results for a num- 
ber of (n,/?)-honeycombs with n G O3 are also reported). 

To summarize, we have studied graphene superlattices 
of defects where a gap at the Fermi level opens because 
of symmetry. Tight-binding and density-functional-theory 
calculations show that the gap is indeed sizable, thereby 
suggesting that the proposed structures may play a role in 
designing future carbon-based electronic devices. 
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